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RECONSTRUCTION OF TRANSIENT ACOUSTIC RADIATION FROM A FINITE 
OBJECT SUBJECT TO ARBITRARILY TIME-DEPENDENT EXCITATION 

BACKGROUND OF THE INVENTION 

[0001] Analyses of transient acoustic radiation are often encountered in engineering 
applications because most vibrating structures are subject to impulsive or transient force 
excitations. In many cases, the transient excitations are unspecified and therefore the structural 
vibration responses are unknown, which makes the analysis of resulting acoustic field difficult. 
While there are available numerous research papers on predicting transient acoustic radiation 
from a finite object subject to the given excitations or time-dependent boundary conditions, 
literature on reconstruction of transient acoustic radiation is almost nonexistent. The closest work 
is the use of spherical near-field scanning to predict far-field acoustic radiation in the time 
domain. This scarceness is caused by a lack of an effective methodology to reconstruct transient 
acoustic radiation based on a limited number of acoustic pressure measurements. 

[0002] One way of dealing with transient acoustic radiation from an arbitrary object is 
via the Kirchhoff-Helmholtz integral theory, which correlates the field acoustic pressure to the 
surface acoustic pressure and normal surface velocity. These surface acoustic quantities are 
governed by the Kirchhoff-Helmholtz integral equation and must be solved before the field 
acoustic pressure can be calculated. For an arbitrary surface, this Kirchhoff-Helmholtz integral 
theory can only be implemented numerically through boundary element method (BEM) by 
discretizing the surface into elements and interpolating the acoustic quantities on the discrete 
nodes. Since the acoustic quantities on these nodes have different emission times for an observer 
in a fixed coordinate system, the integrals become time dependent. Accordingly, one must 
discretize the integrals in both time and spatial domains, thus making numerical computations 
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extremely time consuming. Needless to say, it will be even more difficult, if possible at all, to 
apply the Kirchhoff-Helmholtz integral theory to transient reconstruction. 

[0003] An alternative in transient NAH is to reconstruct the acoustic quantities in the 
frequency domain first, and then take an inverse Fourier transform to retrieve time domain 
signals. Transient acoustic radiation from a vibrating object has been reconstructed in this 
manner. It has been demonstrated that the transient acoustic field can be reconstructed in the 
frequency domain using the Helmholtz equation least squares (HELS) method, and then taking 
an inverse Fourier transform to recover the time history of the normal surface velocity response. 
The infinite integral in the inverse Fourier transform is calculated by direct numerical integration 
for a vibrating sphere in a free field. Needless to say, the numerical computations involved in this 
reconstruction process are huge. Therefore, a direct calculation of an inverse Fourier transform to 
reconstruct transient acoustic radiation is not advisable. 

[0004] The temporal acoustic field has been reconstructed via the so-called non- 
stationary spatial transformation of sound field (NS-STSF). NS-STSF is based on the time 
domain holography (TDH) that processes the acoustic pressures measured by a two-dimensional 
microphone array with the neighboring microphones separated by half wavelength. TDH "can be 
seen as a sequence of snapshots of instantaneous pressure over the array area, the time separation 
between subsequent snapshots being equal to the sampling interval in A/D conversion. Similarly, 
the output of TDH is a time sequence of snapshots of a selected acoustic quantity in a calculation 
plane parallel to the measurement plane." So what one sees are acoustic pressure and intensity 
distributions in the frequency domain at some fixed instants over the measurement time period. 
As such, NS-STSF is not transient NAH in the sense that it cannot provide a three-dimensional 
image of an acoustic signal traveling in space and time. Moreover, NS-STSF is very restrictive: it 
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only allows reconstruction of the acoustic quantities on the same grids as the measurement grids 
on a plane parallel to the measurement plane in a free field. 

SUMMARY OF THE INVENTION 

[0005] In this patent application, we show that transient NAH formulations can be 
derived using HELS method to visualize a 3D image of an acoustic signal traveling in space and 
time based on a finite number of acoustic pressure measurements in the field. It is worth 
mentioning that HELS is valid for arbitrary source surfaces and has been used successfully to 
reconstruct time-harmonic acoustic radiation in both exterior and interior regions. 

[0006] First general formulations for reconstructing transient acoustic quantities on an 
arbitrary surface in a free field are presented. Then general formulations for the temporal kernels 
to be used in the transient NAH formulations are given. Stability of the solutions thus obtained is 
then discussed. Numerical examples of reconstructing transient acoustic radiation from various 
objects are then demonstrated. 

[0007] BRIEF DESCRIPTION OF THE DRAWINGS 
[0008] Other advantages of the present invention will be readily appreciated as the same 

becomes better understood by reference to the following detailed description when considered in 

connection with the accompanying drawings wherein: 

[0009] Figure 1 is a schematic of measurements around an impulsively accelerated 

sphere. 

[0010] Figure 2 is a comparison of the exact and reconstructed normal surface velocity 
distributions on the surface of an explosively expanding sphere. 
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[0011] Figure 3 shows the transient acoustic pressure fields generated by an explosively 
expanding sphere. 

[0012] Figure 4 is a comparison of the exact and reconstructed normal surface velocity 
distributions on the surface of an impulsively accelerated sphere. 

[0013] Figure 5 shows the transient acoustic pressure fields generated by an impulsively 
accelerated sphere. 

[0014] Figure 6 is a comparison of the exact and reconstructed normal surface velocity 
distributions on the surface of an impulsively accelerated baffled sphere. 

[0015] Figure 7 shows the acoustic pressure fields generated by an impulsively 
accelerated baffled sphere. 

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 
[0016] The acoustic pressure reconstruction system 20 of the present invention is shown 
schematically in Figure 1. The system 20 generally comprises a plurality of transducers 22, such 
as microphones, connected to a signal analyzer 24. Alternatively, a digital sound processing 
computer board or software algorithm could be used as signal analyzer 24. The system 20 
further includes a computer 26 receiving and analyzing the data from the signal analyzer 24. As 
is well-known, the computer 26 includes a processor 27 operating a computer programs stored 
on computer storage media 28, which may be one or more of RAM, ROM, hard-drive, CD- 
ROM, DVD, optical or magnetic media, or any other computer-readable medium. Computer 
media 28 includes a computer program, which when executed by processor 27 performs the steps 
described below, including performance of the algorithms of the present invention described 
below. 
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[0017] In Figure 1, nine transducers 22 are placed adjacent to the sound source 42, 
preferably inside an anechoic chamber (not shown). The number of transducers 22 used depends 
upon the complexity of the shape of the object and the accuracy desired, but at least two 
transducers 22 are required. 

[0018] Generally, in operation, the transducers 22 measure the frequency and amplitude 
of sound from the source 42. The gathered data is sent to the signal analyzer 24 which indicates 
amplitude as a function of frequency. This frequency/amplitude data is sent to the computer 26. 
The computer 26 utilizes an inventive method, which will be more fully explained below, to 
determine the transient acoustic pressure fields generated by the source 42. The computer 26 
may display the transient acoustic pressure fields on a display 44. The present invention can be 
used to diagnose general noise sources and transmission since in engineering most structures are 
subject to arbitrarily time-dependent excitations. 

Transient HELS Formulations 

[0019] Consider transient acoustic radiation from a sphere in a free field of density p 0 

and speed of sound c . Assume that the sphere is excited at t = t 0 and prior to that the acoustic 

pressure at any point is zero, i.e., p(x\t) = 0 at t < t 0 . Our objective is to reconstruct the transient 

acoustic quantities everywhere, including the source surface, based on acoustic pressure signals 
measured on a conformal surface T around the source. To this end, let us first define the Laplace 
transform pair for a function /(x;r) as 




o 



and 




(1) 



5 



Docket No. 67021-005 
03-667 



[0020] where s is a complex variable in the Laplace or s -domain, £ is a real constant 
that is greater than the real parts of all singularities of F(x;s). Note that the function /(x;r) in 

oo 

Eq. (1) satisfies the condition j|/(x;f)e~ a |df <<*> for some finite real value of a , and e is 

o 

strategically placed to ensure that the real parts of all singularities of F(x\s) fall on the left-half 
s -plane so the result is bounded as Error! Reference source not found - Accordingly, we can 

express the acoustic pressure in the s -domain as an expansion, which can be written under the 
spherical coordinates as 

N n 



P{r,6,0;s) = X (py m (6,0)C nl (s), (2) 

[0021] where j3-isr/c , h n {0) are the spherical Hankel functions of order n of the first 
kind, Y l n {d,(l>) are the spherical harmonics, and the expansion coefficients C nl {s) are specified by 
solving an over-determined system of equations obtained by matching the assumed-form 
solution (2) to the acoustic pressures measured on a surface enclosing the source at 

x r = {^} T e T , m = 1 to M , where M > J , here J = {1 + N) 2 indicates the maximum number 
of expansion terms in (2). 

P(x r ;^)=^(x r ;.)c(.), (3) 
[0022] where P(x r ;s) and C(s) represent the column vectors containing measured 
acoustic pressures and expansion coefficients, respectively, and H*(x r ;s) is the matrix whose 
elements are given by The expansion coefficients C(s) in 

(3) can be solved as 

C{s)=v(x r ;s}F(x r ; S ), (4) 
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[0023] where <p(x r ; sj = \v(x r ; s f ^(x r ; ¥(x r ; sf is known as pseudo inversion of 
Y(x r ;s) with 

[0024] a superscript H indicating a conjugate transpose. Once C(s) are specified the 
acoustic pressure on a source surface x 0 e S Q or anywhere else can be reconstructed, 

P(x 0 ;^)-G^x 0 |x r ;5)p(x r ;^ (5) 

[0025] where G p (x 0 |x r ;s)= x V(x Q ;s) x lf(x r \sJ is the transfer matrix for reconstructing 
acoustic pressure. 

[0026] The normal surface velocity of a source can then be obtained via the Euler's 
equation 

V(x 0 ;^) = G v (x 0 |x r ;^)p(x r ;^), (6) 
[0027] where G v (x 0 |x r ;5-)= <E>(x 0 ;^)^(x r ; 4 y) t is the transfer matrix for reconstructing the 

normal surface velocity, here Q>(x 0 ;s) = (ip 0 c)~ l d*¥(x 0 ;s)/dt] 0 with r] 0 being the unit normal on 
a source surface and p 0 c the characteristic impedance of the fluid medium. 

[0028] The time-domain acoustic pressure and normal component of velocity on a source 
surface x 0 = (r Q ,0 o ,0 o ) can be obtained by taking an inverse Laplace transform of (5) and (6), 
which is expressible as a convolution in time domain. 

y=i m=l 
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[0029] where jc^eT, m = l to M , and the temporal kernels g; m (* 0 [*I;') and 
Sm^I**)^ given by 

g^t|^^)=^ Jg;(x 0 ^;.)^, (9a) 



[0030] Needless to say, the temporal acoustic quantities anywhere in the field can be 
obtained by substituting the coordinates at x 0 in Eqs. (7) and (8) by those at x = (r,0,^) . 
Evaluation of Temporal Kernals 

[0031] The integrands in Eqs. (9a) and (9b) are given by 
^(^K;^^!^^)^^^;^ and G;Jx 0 ^;.)=O y .(x 0 ;.)y.(^;^ 5 and the functions 

Y y (r f 5) S ^(r,^^j) = ^09)i;'(ftrt and S> , ■* - (r i j) = (^c)" l 3/a7kW^(ftrt]. where 

the indices j , n , and / are related together through j = n 2 + n + 1 + 1 , with n increasing from 0 

to N and / from -n to/2. 

[0032] To facilitate evaluations of the temporal kernels, we recast the infinite integrals in 
Eq. (9) into contour integrals. 

g%l^lt)«±jGl(^;s)erds, (10a) 
^UK^)=-^^(^;.) e ^, (10b) 



t 



[0033] where the contour £ consists of a straight line from e-iR to £ + iR and a 
semicircle of radius R in the left-half s -plane in a counterclockwise direction. The contribution 
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of this integral from semicircle is zero as /? — » <*> . This is obvious if we notice that the integrands 
in Eq. (10) contain the factor e -( ct + rr - r ») R ™Yl c that decays exponentially as R -> oo since t > 0 and 
r r > r Q , where r r and r 0 are the radii of source S 0 and measurement F surfaces, respectively, 
and y is the polar angle with respect to the real axis for a point on the semicircle. The contour 
integrals in (10) can be evaluated using the residue theory. To this end, we replace the spherical 
Hankel functions and their derivatives in G^(jc 0 |jc£;s) and GJ m (x 0 |*^;.s) by polynomials of s 



[n + q)\ 



<7=0 1 



2«q\{n-q)\r) 



j_ 

,4+1 



(11a) 



h (ir \ _ ,-„/ c f (n + qWrs/c)+(q+l)] (cT l 1 (n ) 

k {r)- e U r*4(n- q \ lr) *«• (Ua) 

[0034] Substitute Eq. (11) into G^(x 0 |jc£;s) and GJ m (* 0 |x^;s), use partial expansion to 

recast the polynomial into simple poles, and sum them up in terms of simple real and complex 
poles. Doing so, we obtain 



G p Jxfas)= Vfa.sfrfas) = e^< 



G] m {xfas)= * fas)* fas) = J*-** 



+ e {£-r 0 )s/c 



fit- 1 



, (12a) 



£o h s 2 -2Rc{s%) s + \s; 



2 
<? J 



*l0 *r 



9=0 S 



, s + \s v : e 



(12b) 
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[0035] where (sjjj^ and (sj m ) q denote the qth simple real and complex poles for 

reconstructing surface acoustic pressure, respectively, n p jm (s^ m )^ and tl p m {s^ m )^ stand for the 

residues that correspond to (s^) and (sj m ) q , respectively, Q p jm and Q p jm are total numbers of 

simple real and complex poles for Gj^Jt 0 |jc£;s). Similarly, (s]%) and {s v ^) q indicate the qih 
simple real and complex poles for reconstructing the r -component of the surface velocity, 
respectively, 11]^ (jj^) and n^(s^)^ represent the residues that correspond to and 

[s V jm) q » respectively, and and are total numbers of simple poles for the r -component of 
GJ m (jc 0 |x^;^). The gth real and complex poles, residues, and numbers of poles for reconstructing 
the 0 -component of the surface velocity are defined as (s v *) q9 ($Jf) fl> n v £(s v £) g , A-ffe 
Q V j£ , and Q£ , respectively. The q th real and complex poles, residues, and numbers of poles for 
reconstructing the 0 -component of the surface velocity are given by (sjjf) , {s]m) q > n^(j^)^, 
njfel , Q'J , and Q v £ , respectively. 



[0036] The general solutions for the temporal kernels can now be obtained by 
multiplying sums of the residues by iln . 
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(13b) 



[0037] where H{t - r T m - rjc) is the Heaviside step function. 

[0038] Equation (13) illustrates that the temporal kernels consist of the components of 
weighting functions. In particular, the first group of terms in the square brackets of Eq. (13) that 
result from simple real poles (sj m ) q , [ s ]m) q > { s 7m) q > anc * { s 7m\ decay exponentially in time. 

While the second group of terms resulting from simple complex poles (sj m ) q > {^]m\ > \ > anc * 

[sfa) oscillate sinusoidally, with their amplitudes decaying exponentially in time and decreasing 

as the angular frequency increases. The real parts of complex poles Re(^ m )^, Re(^ r )^, 

Rzfyn)^ , and Re^)^)^ dictate the rates at which the weighting functions decay exponentially in 

time, and the imaginary parts Im(^ m )^, Im(sj m r )^, Im(^)^, and lm(s^) depict the angular 

frequencies of the weighting functions. When the temporal kernels given by Eq. (13) are 
convoluted with measured acoustic pressure signals, we see transient acoustic pressure waves 
whose amplitudes decay exponentially in time as they travel in space. 
Stability of Solutions 

[0039] As in all inverse problems, we must face the ill-posedness difficulty in 
reconstructing the temporal acoustic quantities. In time-harmonic NAH, this ill-posedness 
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problem is dealt with via regularization that smoothes the dependence of solution on the input 
data. A common approach is to employ Tikhonov regularization (TR) with certain parameter 
choice method together with singular value decomposition (SVD). The objective of TR is to 
strike a balance between minimal residues and minimal solution norms. The purpose of SVD is 
to identify and eliminate the small singular values that may amplify the errors embedded in the 
input data and cause distortion in the reconstructed images. 

[0040] In transient NAH regularization can be accomplished by eliminating the simple 
poles that stay close to the imaginary axis. Since the contour integral is carried out in the left-half 
s -plane, the real part of a simple pole is negative and the amplitude of the corresponding 
temporal kernel decays exponentially in time. This guarantees the stability of temporal solutions. 
However, when a simple pole approaches the imaginary axis, its real part diminishes and the 
amplitude of the corresponding temporal kernel decays very slowly. Under this condition, the 
errors embedded in the input data remain finite and get transferred to the reconstructed acoustic 
quantities. 

[0041] To circumvent this difficulty, we propose a low-pass filter to eliminates all simple 
poles whose real parts are larger than certain tolerance € . In this way, the amplitude of the 
acoustic pressure remains finite as the wave propagates in space and time. Accordingly, we can 
rewrite Eq. (12) as 
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h ^ST + ^, 2 -2Re(^) /+ |^ 



9 ' 



(14b) 



<7=0 



[0042] where © p jm and B p jm are total numbers of the principal simple poles for 



jm 



reconstructing surface acoustic pressure, Q v jm and Q v jm are those of the principal simple poles for 

reconstructing normal surface velocity in the time domain, and the rest simple poles whose real 
parts are larger than e are eliminated. Thus this low-pass filter expands the acoustic field in 
terms of the acoustic modes that are bounded in the time domain, which is equivalent to what 
truncated SVD does for time-harmonic NAH applications. The regularized temporal kernels can 
then be written as 



rm(*o|*m»*)~ 



q=0 



^ Im jf J 1 q 

V j™ In 



H 



f~k— k ,(15a) 



5 



(*oK;') = 



e v _-i 



9=0 



6 M*U 



( r - "\ 



(15b) 



Numerical Examples 

[0043] In this section, we use the transient NAH formulations (7) and (8) derived in 
above to reconstruct the transient acoustic quantities generated by a sphere. To facilitate 
validations of the reconstructed field, we consider a sphere of radius r 0 = 1 m subject to an 



impulsive excitation for which the analytic solution is readily available. 
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v(r 0 ,fl 0 ^;') =<! fc^o4W'-0. ( 16 ) 

[0044] where v(r o ,0 o ,0 o ) is the unknown amplitude of the normal velocity on the source 

surface. 

[0045] In these examples, the measured acoustic pressures are generated analytically 
along a line segment that intersects at 135° with respect to the a: -axis with a minimum distance 
of d = 0.1 m (see Figure 1) These acoustic pressures are taken as the input to Eqs. (7) and (8) to 
reconstruct the surface acoustic pressures and the normal surface velocities in the time domain, 
which are then compared with the analytic results. 

Explosively expanding sphere 

[0046] ^ e fj rst exam pj e d ea i s with reconstruction of the transient acoustic field 
generated by an explosively expanding sphere in free space. One can view such a transient 
acoustic field as being produced by an explosion. Since the field is spherically symmetric, the 
dimensionless acoustic pressure depends only on the radial distance and time and the analytic 
solution is given by 



gM _ 



ae 



~(ct-r)/a 



-H 



t — 



(17) 



p 0 cv r \ c) 

[0047] Equation (17) is used to generate 24 field acoustic pressures, which are taken as 
input to Eqs. (7) and (8) to reconstruct the transient acoustic field. Test results show that 
solutions given by Eqs. (7) and (8) converge very rapidly. This is no surprise since the acoustic 
field is produced by the zeroth order of the spherical Hankel function. As a result, satisfactory 
reconstruction can be obtained without using regularization. 

[0048] Table 1 illustrates comparison of the first eleven reconstructed expansion 
coefficients and the exact ones. For consistency, we show comparison of first eleven expansion 
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coefficients in all cases. Results show that the transient acoustic field produced by an explosively 
expanding sphere in a free field can be reconstructed using one expansion function (or taking one 
measurement). 



r 
*j 


Reconstructed values 


Analytic Values 


Ci 


9.99999E-01 


1.0 




2.16846E-16 


0.0 


£ 3 


5.99974E-16 


0.0 


c 


-2.41091E-16 


0.0 


Cs 


3.5861 1E-16 


0.0 


Ce 


2.62832E-16 


0.0 




-5.34793E-16 


0.0 


C s 


-4.50502E-16 


0.0 


Co 


4.49696E-16 


0.0 


Co 


1.87202E-16 


0.0 


Cn 


-3.71210E-16 


0.0 



[0049] Table 1. Comparison of coefficients ^ for an explosively expanding sphere. 

[0050] Figure 2 exhibits the comparison of the reconstructed and exact normal surface 
velocities. Once the normal surface velocity is determined, the radiated acoustic pressure at any 
instant and any location can be reconstructed. Figure 3 shows three-dimensional images of the 
instantaneous acoustic pressure distribution at t = 3.24 (ms), 4.41 (ms), 5.88 (ms), and 7.35 (ms) 
resulting from a sphere that explodes at t = 0 (sec). 
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Impulsively accelerated sphere 

[0051] The second example concerns reconstruction of the transient acoustic field 
resulting from an impulsively accelerated sphere in a free field. One can view this as a crack on a 
baseball. The analytic solution for the temporal acoustic pressure generated by this sphere is 
given by 



p(r>0\t) _ flcosfl c -{ct-r)ia H 
p 0 cv r 



( A 
t — 

V cj 



cos 



ct-r\ I a 



r) 



sin 



ct-r 



(18) 



[0052] where co$0 = ti*e R , here is the unit vector in the direction of wave 
propagation from the source to the receiver. As in the first case, we collect 24 input data and use 
them in Eqs. (7) and (8) to reconstruct the acoustic field. Table 2 shows comparison of the 
reconstructed expansion coefficients and the exact ones for an impulsively accelerating sphere. It 
is seen that in this case all coefficients are zero except the second one and excellent agreement is 
obtained. 





Reconstructed values 


Analytic Values 


Ct 


3.98736E-17 


0.0 


<r 2 


9.99999E-01 


1.0 


c, 


-1.06066E-16 


0.0 


c 


3.24818E-11 


0.0 


Cs 


1.97787E-16 


0.0 


Co 


-3.32822E-10 


0.0 


Ci 


-1.05409E-17 


0.0 




3.70379E-10 


0.0 
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-1.60172E-16 


0.0 


<r I0 


-1.02017E-11 


0.0 




8.92684E-17 


0.0 



[0053] Table 2. Comparison of coefficients Cj f° r an impulsively accelerating sphere. 

[0054] The reconstructed normal component of velocity is displayed in Figure 4, which 
agrees very well with the analytic value. Numerical results indicate that theoretically we can 
reconstruct the transient acoustic field generated by an impulsively accelerated sphere using four 
expansion functions only, which correspond to four measurements. This is because the first 
expansion term represents a monopole and the next three depict dipoles in (x,;y,z)-axis 
directions, respectively. Figure 5 illustrate the instantaneous acoustic pressure wave fronts at / = 
3.24 (ms), 4.41 (ms), 5.88 (ms), and 7.35 (ms) from the sphere impulsively accelerated at t = 0 
(sec). 

Impulsively accelerated baffled sphere 

[0055] In this example, only a small portion of the surface is impulsively accelerated and 
the rest area is rigidly baffled. The resulting acoustic pressure field can be viewed as being 
produced by popping a champagne bottle in a free field. Mathematically, the normal component 
of the surface velocity is expressible as, 

v(r 0 ,^0 = ^o(^o^)[^(^ + ^o)-^(^-^oM--^o) > (19) 
[0056] where 20 o is the vertex angle (# 0 =30°) and v(r s ,0) is to be reconstructed. 

[0057] The analytic solution for the radiated transient acoustic pressure is given in terms 
of an infinite series. The same procedures as described above are followed in reconstruction. 
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[0058] Table 3 displays comparison of analytic expansion coefficients and reconstructed 
values. Because of a rich content of high spatial frequencies, which is reflected in the presence of 
a sharp edge in the normal surface velocity given by Eq. (19), a large number of expansion 
functions are needed to approximate the velocity field. Figure 6 displays the comparison of the 
exact normal surface velocity and reconstructed ones using different number of expansion 
functions. It is seen that solutions given by Eq. (8) converges to the exact normal surface velocity 
and the accuracy of reconstruction increases with the number of expansion functions. 
Theoretically, if input data are complete and error free, the reconstructed normal surface velocity 
may converge to the true value as the number of expansion functions approaches infinity. 
However, this will never happen because in reality the measured data are always incomplete and 
inaccurate. Therefore, we rely on regularization to truncate the matrix in Eq. (8) so that errors in 
input data do not get amplified in the reconstruction process. 





Reconstructed values 


Analytic Values 


& 


6.69999E-02 


6.70000E-02 




1.87500E-01 


1.87500E-01 




2.70633E-01 


2.70633E-01 




3.00783E-01 


3.00781E-01 


Cs 


2.74016E-01 


2.740 16E-01 


£ 6 


1.98730E-01 


1.98730E-01 


<T 7 


9.3453 1E-02 


9.34529E-02 


<Ts 


-1.76239E-02 


-1.76239E-02 


to 


-1.10301E-01 


-1.10301E-01 
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-1.65868E-01 


-1.65869E-01 




-1.75139E-01 


-1.75140E-01 



Table 3. Comparison of coefficients ( } for a partially impulsively accelerating baffled sphere. 

[0059] Once the normal velocity on the surface is determined, the entire acoustic 
pressure field can be reconstructed using Eq. (7). Figure 7 shows the transient acoustic pressures 
produced by an impulsively accelerating baffled sphere at different instants. Results indicate that 
the pressure pulse is produced primarily around the area where a velocity pulse is initiated and 
the amplitude becomes extremely low beyond certain polar angle. 

Conclusions 

[0060] General algorithms are developed for reconstructing transient acoustic radiation 
from a finite object in a free field. While examples shown are for a sphere, these algorithms are 
valid for an arbitrarily shaped object subject to arbitrarily time-dependent excitations. This is 
because the expansion functions used in these algorithms to describe the acoustic field are 
complete. These algorithms allow for an efficient and accurate reconstruction of a transient 
acoustic signal as it travels in space and time. 
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